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Long range order in non-equilibrium interacting quantum spin chains 
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We conjecture that non-equilibrium boundary conditions generically trigger long range order in 
non-equilibrium steady states of locally interacting quantum chains. Our result is based on large 
scale density matrix renormalization group simulations of several models of quantum spin 1/2 chains 
which are driven far from equilibrium by coupling to a pair of unequal Lindblad reservoirs attached 
locally to the ends of the chain. In particular, we find a phase transition from exponentially decaying 
to long range spin-spin correlations in integrable Heisenberg XXZ chain by changing the anisotropy 
parameter. Long range order also typically emerges after breaking the integrability of the model. 
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Non-equilibrium physics of low dimensional locally in- 
teracting quantum many-body systems is becoming an 
increasingly hot topic, in particular due to enormous re- 
cent progress in experimental techniques with ultra-cold 
atomic gases and optical lattices [l| , and arguable connec- 
tion to important open challenges in condensed matter 
physics such as high-temperature superconductivity and 
solid-state quantum computation. While thermalization 
in isolated many-body systems seems to be well under- 
stood 13 , and near equilibrium physics is essentially given 
by linear response formalism, far from equilibrium even 
some very basic concepts remain unclear. 

One of the key questions that we address in this Let- 
ter is: Driving a large one-dimensional locally interacting 
quantum system out of equilibrium by coupling to two 
unequal reservoirs at the system's ends, will the non- 
equilibrium steady state (NESS) in the bulk be well de- 
scribed by the notion of local equilibrium? If yes, then the 
properties of NESS are essentially given by local Gibbs 
state with a well defined local temperature (and perhaps 
local chemical potential) field(s). If not, then the local 
temperature has no meaning, and the physics of far from 
equilibrium system is qualitatively different than in equi- 
librium. While Mermin-Wagner's theorem [3] prohibits 
certain types of order at finite-temperature equilibrium, 
(arbitrarily weak) non-equilibrium driving at the bound- 
ary of the system breaks not only the translational invari- 
ance but typically also the symmetries of the Hamiltonian 
and opens the possibility of long range order (LRO). As 
has been shown for an open XY spin 1/2 chain, LRO 
may set in abruptly by varying the system's parameters 
through a non- equilibrium quantum phase transition [j|. 

Nevertheless, since the XY chain is equivalent to a 
quasi-free fermionic model, one might argue that LRO 
might be a non-generic effect which would typically not 
survive switching on the interaction, or even more so, 
breaking the integrability of the model. In this Letter we 
show that this is not the case by performing large scale 
density-matrix-renormalization-group (DMRG) simula- 
tions of NESS for integrable and non-integrable versions 



of anisotropic locally interacting Heisenberg (XXZ) spin 
1/2 chain. By changing the anisotropy parameter we find 
a clear evidence of non-equilibrium quantum phase tran- 
sition from exponentially decaying spin-spin correlations 
to LRO where spin-spin correlations saturate and do not 
depend on the system size. Perhaps even more strik- 
ingly, we find that LRO is typical for non-integrable de- 
formations of the model. Note that the physics of NESS 
can be very different from the much studied quench dy- 
namics where one can observe at most a transient non- 
equilibrium phenomena. 

The model we study is an anisotropic nearest-neighbor 
interacting spin 1/2 Heisenberg XXZ chain 
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The Hamiltonian ([T]) is an important model for several 
reasons. It describes various spin-chain materials [5] like 
cuprates and as a particularly promising experimental 
technique, it can be realized with cold atoms in an opti- 
cal lattice [6|. From the theoretical side, ([T]) is perhaps 
the simplest model of strongly interacting electrons (to 
which it can be mapped by Jordan- Wigner transforma- 
tion). Understanding the Heisenberg model would there- 
fore provide a guidance to physics of more realistic mod- 
els of interacting electrons like e.g. the Hubbard model. 
Despite integrability of the Heisenberg model J7| , cal- 
culating quantities beyond, for instance, the ground state 
energy or low-energy excitations is exceedingly difficult. 
Spin-spin correlation functions are particularly impor- 
tant quantities. Regarding equilibrium correlation func- 
tions in the Heisenberg model the following is known 
7]: (i) at infinite temperature there are no correlations, 
(ii) at finite temperature correlations asymptotically de- 
cay exponentially with distance, (iii) at zero tempera- 
ture (ground state) correlations decay asymptotically as a 
power-law of the distance in the gapless phase for | A| < 1, 
while they decay exponentially in the gapped phase of 
|A| > 1. For non-equilibrium stationary states almost 



nothing is known about the correlation functions, some 
preUminary hints can be found in [8|. The aim of this 
work is to focus on long range correlations in NESS. 

While the properties in close-to-equilibrium situation 
can be studied using a linear-response formalism, far from 
equilibrium one is forced to consider genuine nonequilib- 
rium setting where the central system is coupled to reser- 
voirs at different potentials, for instance, at the chain 
ends. We are going to take into account the Markovian 
reservoirs in an effective way by using the Lindblad equa- 
tion. Time evolution of systems's density matrix is there- 
fore governed by the Lindblad master equation [9| 
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The dissipator £!°^^^ — £}^ ^ £,^ is a sum of two parts, 
one acting on the left end and the other acting on 
the right end of the chain. Each one is written in 
terms of a set of Lindblad operators L^ as C^{p) = 
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L, R. Most of the 



time we shall use the one-spin bath for which Lindblad 
operators act only on a single spin. Two Lindblad opera- 



tors at each end are i^2 = \/l ^±^1 ^ ^? 
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where a/, = h{(T^ ± ia^). In numeri- 



l±tanh£f^ 

cal simulations we always consider strong coupling 7=1. 
These operators are such the the stationary state of the 
(say, left) bath spin only, C^{pi) = 0, is pi oc exp (fJ-hcrl) 
with polarization (cr^) — tanh/iL. One can argue that 
such a situation formally corresponds to an infinite tem- 
perature with fixed bath magnetizations. Later we shall 
also consider simple two-spin baths where two border 
spins are coupled to a bath in order to simulate the effect 
of finite temperature. The choice of the bath is qualita- 
tively inessential for the results shown. 

To find NESS (fixed point of the Liouvillean) we sim- 
ulate the Lindblad equation ^ using time-dependent 
DMRG (tDMRG) [iBl with matrix product operator 
(MPO) ansatz. This enables us to study systems much 
larger than it would be possible with brute-force diago- 
nalization. After sufficiently long time, the density ma- 
trix p{t) converges to a time-independent NESS. Once 
NESS density matrix pness is obtained in terms of MPO, 
various expectations can be efficiently evaluated. [l8| 

We now turn to the study of the main object of the 
present work, the 2-point spin-spin correlation function 

C{i,j)^ trpNESSo-jO-j - (tr/9NESsCTD(trpNESS0-j)- (3) 

C{i,j) is usually found to depend on the relative dis- 
tance r = |i ~ j| ranch stronger than on the perpendic- 
ular coordinate, thus we often show a one-dimensional 
cross-section of C{i,j) along the skew-diagonal, that is 
C(r) := C(fc, k + r), where k = [s±ij - [§J . Sometimes, 
when larger fluctuations among the values of C{i,j) are 
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FIG. 1: Spin-spin correlation function C{r) in open XXZ 
model at A = 0.5 (a) and A — 1.5 (b), both for fiL.n = 
±0.22. Symbols are from exact calculation, dashed curves 
from tDMRG, while full curves suggest n — >■ 00 asymptotics. 



noticed, we average the absolute value of the correla- 
tion over all pairs of points with a constant distance 
(|C(r)|) := {\C{k, k -\- r)\)k, and in order to characterize 
LRO, as order parameter we compute the residual cor- 
relator as an aditional average over r between n/4 and 
3n/4, (|q):=(|C(fc,fc-Hr)|)^/^^'■=^^"/^ 

It turns out that the spin-spin correlation function in 
NESS depends strongly on the anisotropy A. For exam- 
ple, for A = 0.5 the correlation function (its cross sec- 
tion C{r) is shown in Fig. [1^) shows exponential decay, 
with asymptotic scaling C{r) ^ e"''/^, with ^ = 2.0. We 
can see that the correlation decay rate ^ is independent 
of the chain size. Therefore, correlations in NESS are 
short-ranged. However, after increasing the anisotropy 
the situation becomes dramatically different, for exam- 
ple, for A = 1.5 the correlator C{r) clearly exhibits LRO 
as shown in Fig. [T]d. With increasing the chain size n 
the correlation function converges to the limiting form 
for which we can observe several things: (i) spatial de- 
pendence is a function of the scaled variable r/n only, 
for large n, meaning that there is LRO in NESS, (ii) 
the size of the correlations scales with the driving as 
^ (Ap)^ = (^L — Mr)^ and is independent of n. Further- 
more, correlations are positive (the point where they be- 
come negative moves towards r/n — ;■ 1 as n is increased). 
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FIG. 2: Residual correlator (|C|) in open XXZ model against 
the scaling variable (A — Ac)n with Ac = 0.91, for n = 20 
(circles), 30 (diamonds), 40 (squares). The inset shows un- 
sealed data. Here /^l.r = ±0.22 whereas practically the same 
scaling of {\C\)/{Afi)'^ was obtained also for /il.r = ±0.05. 



The limiting form of the correlation function is perhaps 
close to C{i = nx,j = ny) ~ (A/i)^a;(l — y). Very similar 
form of the 2-point correlation function, however with a 
size dependent prefactor (l/n), has been found in cer- 
tain non-interacting, dissipative (noisy) or even classical 
models: in a ballistic exactly- solvable XY model 4,], nu- 
merically observed in the diffusive Heisenberg model with 
dephasing 11] and analytically explained in a diffusive 
XX model with dephasing [12| , but also in some classical 
lattice gasses like simple exclusion processes 13| or its 
quantum analogs 141]. The fact that long range corre- 
lator here does not decay with n perhaps hints at truly 
quantum coherent character of these correlations. 

From Fig IT] we can conclude that somewhere between 
A = 0.5 and A ~ 1.5 the range of correlations must 
change from short to long range, i.e., there is possibility 
of a nonequilibrium phase transition. To find if this is 
the case and where is the transition point we study the 
residual correlator (|C|) against parameter A for several 
system sizes. Remarkably, numerics (shown in Fig. [2]) 
indicates the scaling (|C|) ^ /((A — Ac)n''), where / is 
some universal scaling function and t/ « 1.0, suggesting 
that the transition becomes non-continuous in the ther- 
modynamic limit. Results seem to be accurate enough to 
support the claim that the observed critical anisotropy 
Ac « 0.91 is distinct from the zero temperature critical 
anisotropy A* = 1 where the ground state gap opens. 

One may argue that LRO observed above might be due 
to integrability of the XXZ model. In order to check the 
effect of (non)integrability we add a staggered magnetic 
field to the Hamiltonian Q, namely H' = H+X]?=i ^j"'? 
with h2j = h, h2j+i = 0, for h = —0.5, which is known to 
bring the model in the regime of quantum chaos [8| . Us- 
ing the same driving as before, we find that NESS again 
exhibits long range correlations. In order to somewhat 
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FIG. 3: Smoothed spin-spin correlation function (|C(r)|) in 
the non-integrable deformed XXZ model with A — 0.5 and 
staggered transvere magnetic field h = —0.5, /iL,R = ±0.1. 



isolate the effect of integrability breaking alone, we take 
A = 0.5, which in the integrable limit exhibited short 
range correlations. In order to smoothen the oscilla- 
tions in the correlation function due to staggered field we 
plot in Fig. [3] the average correlator (|C(r)|) and clearly 
observe LRO which is again consistent with a constant 
residual correlator (|C|) — 0{n^). 

All results presented so far have been for symmetric 
driving, fj,^ = — /zl, which results in NESS having aver- 
age magnetization zero and the average energy density 
also zero. This can be interpreted as NESS at an infi- 
nite average temperature, although care must be taken 
when discussing the temperature in an integrable open 
system because thermalization is absent for un-adjusted 
reservoirs due to lack of ergodicity [l5| . Thus the struc- 
ture of the ground state very likely plays no role in the 
above NESS correlations. An interesting question is how 
do the correlations in NESS change as we lower the ef- 
fective temperature? Also, correlations could depend on 
the average magnetization, as does for instance the spin 
transport in a linear response regime. In the zero mag- 
netization sector there appears to be a transition from 
ballistic spin transport for |A| < 1 to a diffusive trans- 
port in the gapped regime |A| > 1. On the other hand, 
for nonzero magnetization spin transport is always bal- 
listic as exemplified by a nonzero Drude weight 16|. To 
study correlation function in NESS at a nonzero energy 
density (i.e., at finite "temperature") we used C^'^^^ cou- 
pling to two spins at the boundary or a one-spin bath 
with non-symmetric driving. Similar results are obtained 
in both cases and we are going to show only the ones 
for a two-spin bath as the relaxation to NESS there is 
faster. Lindblad operators on two boundary spins are 
chosen in such way that the invariant state of these two 
spins would be a grand-canonical state at some target 
temperature and inagnetization in the absence of the 
bulk Hamiltonian [ll|, |l7|. Effective temperature field 
of NESS is then prescribed by equating the measured 
local magnetization and energy density in the bulk to 




FIG. 4: Finite temperature correlation function C{r) for open 
XXZ model at A = 1.5. Frame (a) is for T^eas « 33 while (b) 
is for Tmoas ~ 5.8 (both measured in the middle of the chain). 
At short distances (r < 10 in (b)) we have an exponential 
decay of C{r) with the size-independent decay rate. For large 
r, we find the residual correlator scaling as (|C|) ~ (Ap) /n" 
(insets), with a ~ 0.4 in (a) and a « 2 in (b). 



the grand-canonical one at some "measured" tempera- 
ture Tnicas and potential /imoas through canonical den- 
sity matrix p ~ exp [—{H — fj,S^)/T], see [15| for details. 
The results for correlation function at two different mea- 
sured temperatures in the center of the chain are shown 
in Fig.Ufor the gapped XXZ model A = 1.5 at zero stag- 
gered field h = 0. These two results for finite "tempera- 
ture" show that the decay of the correlation function has 
two regimes. For small distances r, we have a short-range 
exponential decay, with the decay rate depending only 
on the temperature and not the chain length. This can 
be interpreted as a finite-temperature equilibrium part of 
the correlation function. On the other hand, at larger 
distances long-range correlations still persist, but scaling 
algebraically with the system size as (|C|) ^ (A/i)^/n", 
with some nonzero power a. This second regime can be 
interpreted as being genuinely nonequilibrium. Noting 
that at an infinite temperature. Fig. [Hd, we have a = 
as discussed before, while at lower temperatures, Fig. |4l 
a increases. We conjecture that a goes from at an 
infinite temperature to a — >■ oo at zero temperature. In- 
creasing exponent a with decreasing temperature might 
be in line with the exponential decay of correlations for a 



gapped ground state in equilibrium. Note however, that 
in the equal-bath (equilibrium-like) situation, A/i — >■ 0, 
all correlations in NESS decay exponentially. 

In conclusion, we have presented strong numerical sup- 
port for the existence of long range correlations in bound- 
ary driven non-equilibrium steady states of locally inter- 
acting quantum spin chains. We have shown the exis- 
tence of quantum phase transition in an open XXZ spin 
1/2 chain with spontaneous emergence of long range cor- 
relations with increasing the anisotropy, and/or breaking 
the integrability of the model. These long range cor- 
relations are particularly strong, as they - evaluated at 
the fixed scaled distance r/n - do not depend in the sys- 
tem size n, and are thus fundamentally different from the 
'hydrodynamic' long range correlations observed for in- 
stance in certain classical far from equilibrium stochastic 
processes [13j . This phenomenon could have far reach- 
ing consequences for far from equilibrium physics in low 
dimensional quantum many-body systems, such as com- 
plete absence of local equilibrium and non-existence of 
local thermodynamic potentials. The work is supported 
by grants Pl-0044, Jl-2208 of Slovenian Research agency. 
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